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ABSTRACT 

This  investigation  was  undertaken  to  determine  the  feasibility  of 
computing  the  friction  velocity  from  available  micrometeorological  data. 
This  could  be  of  great  assistance  in  finding  the  wind  velocity  profile  by 
Deacon's  equation  and  the  eddy  shearing  stress  in  the  lower  portion  of  the 
friction  layer. 

This  investigation  was  carried  out  at  the  U.  S.  Naval  Postgraduate 
School,  Monterey,  California,  during  the  period  February  1959  -  May  1959, 
in  partial  fulfillment  of  the  requirements  for  the  degree  of  Master  of 
Science  in  Meteorology. 
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1.   Introduction 

The  friction  velocity,  a  stability  parameter,  and  a  term  related 
to  the  Richardson  number  were?  computed  using  micrometeorological  data. 
The  equations  used  for  the  computation  of  these  three  parameters  are 
derived.   These  equations  are  unique  in  that  they  combine  the  different 
theories  of  two  well-known  investigators  of  surface  layer  turbulence  into 
an  equation  for  the  friction  velocity  and  an  equation  for  the  wind  profile 
in  the  surface  layer. 

The  micrometeorological  data  used  in  this  investigation  Is  discussed, 
as  well  as  the  instruments  used  In  obtaining  the  data  and  the  observa- 
tional errors.  A  theory  Is  presented  on  errors  of  aerated  temperatures 
close  to  the  ground. 

After  modifying  the  derived  equations  into  a  form  suitable  for  using 
with  an  electronic  digital  computor,  the  so-called  "constants"  in  these 
equations  are  discussed.   The  advantages  and  disadvantages  of  using  an 
electronic  computor  for  this  type  of  an  investigation  are  then  summarized. 

Various  figures  and  graphs  are  used  to  illustrate  the  diurnal  varia- 
tion as  well  as  the  possible  variation  with  height  as  determined  from 
this  Investigation.   A  comparison  of  these  results  with  those  of  previous 
investigators  is  made,  Indicating  close  agreement  in  most  eases. 


2.   Derivations 

Deacon  [2],  in  the  derivation  of  his  wind  profile  equation,  as  well 
as  most  other  well-known  investigators  in  the  field  of  surface  layer 
turbulence,  considered  only  two  virtual  turbulence  parameters  -  namely,  the 
mixing  length  jl   and  the  friction  velocity  u#,  where 


»*-# 


Here  "■-  is  the  eddy  stress  and  a  is  the  air  density.   Lettau  [10]  introduces 

a  third  parameter,  the  vertical  mixing  velocity  w,  that  is  quite  important 

if  nonadlabatic  rather  than  neutral  conditions  are  investigated.  He  also 
has  a  term  which  he  calls  the  friction  velocity  u,  where 


UW  m    -X 
/° 

However,  when  the  friction  velocity  is  mentioned  in  this  paper  it  will  refer 
to  the  first  definition. 

The  surface  layer  is  defined  as  that  layer,  immediately  above  the 
earth's  surface,  in  which  the  eddy  stress  is  constant  with  elevation  through- 
out the  layer.   Eqs.  (1-9)  are  descriptive  of  this  layer.  The  notation  is 
that  employed  by  Lettau  in  [10]. 

In  a  neutral  surface  layer  the  equations  below  are  applicable. 

•      * 
a      a  (1) 


Throughout  this  paper  u  will  be  understood  to  mean  the  wind  velocity 
at  some  level  averaged  over  some  Increment  of  time  -  say  six  minutes. 


Jl        »  k(z  +  z  )     k  »  von  Karman's  constant  (2) 

z  «  roughness  parameter 


* 
a 
a 


4-^*-  (3) 


u  w   s   ^l  (4) 


V? 


• 


»*         1/-^* 


It  should  be  noted  that  in  Eqs.  (1-5),  u  ,  w  ,  and  C  are  constant 

em      8  A 

with  elevation  in  the  neutral  surface  layer. 

In  a  non-neutral  surface  layer  the  next  four  equations  are  used. 

u   *   w  (6) 


uw  m    2L  (7) 

P 

»*  -  if  (8) 

r^  *  r  (9) 

f   .  *      * 

In  [10  J,  Lettau  derives  a  relationship  between  w  and  w  which  is  of 
fundamental  importance  in  the  development  of  the  remainder  of  this  paper. 
Two  basic  assumptions  are  necessary  in  the  derivation  of  Lettau *s  relation- 
ship: 

1.   The  ratio  w/^  does  not  depend  upon  variations  of  thermal  stratifi- 
cation, so  that 


w  •  wa  (10) 


Eq.  (10)  indicates  that  the  vertical  mixing  length  is  changed  in  exactly 
the  same  ratio  as  the  vertical  mixing  velocity  if  the  thermal  stratifi- 
cation varies. 

2.   The  second  assumption  is  that  the  buoyant  acceleration  influences 
only  the  vertical  parameters  of  turbulence.   This  assumption  is  necessary 
in  order  to  fix  absolutely  the  variation  of  mixing  velocity  and  mixing 
length  with  thermal  stratification.   In  an  adiabatic  atmosphere  the  term 
ivwCA   can  be  considered  to  represent  the  vertical  component  of  a  turbulent 
acceleration.   Thus  we  have  the  vertical  acceleration  equation: 

W2  -  Waf  -  J7^^  (11) 

where  g  *  acceleration  due  to  gravity,  dB/dZ    *  change  in  potential  tempera- 
ture with  respect  to  height,  and  T   «  mean  absolute  temperature  averaged 
with  respect  to  height  in  the  layer  considered.   Eq.  (11)  indicates  that 
the  mean  vertical  turbulent  acceleration  in  a  non-adiabatic  surface  layer 
equals  the  adiabatic  acceleration  plus  the  buoyant  acceleration. 
Dpon  multiplying  Eq.  (11)  by  X/vj} 


*  2    A  Tn  v£* 


but   from  Eq.    (10),    Wywt  =  xV{2    ,    so  upon  substituting 


jr  _   x        ^Tz* 


n     a       Tm  3,2. 


multiplying  through  by  ^  /j£ 2        gives 


or 


i  -'  A  _    rn** 


^4 

Let  x  =    az    *  (12) 


T   *2 


Therefore,     x  =  J±  —  /    ,  (13) 


Thus  we  finish  with  Lettau's  derivation  by  putting  Jl/Jl*     in  terms  of  x, 


A   =  _ L_  (14) 


The  parameter  x  is  important  in  the  computations  to  be  made  in  this  paper. 

The  next  derivation  was  originally  developed  by  Martin  [12]  combining 
Deacon's  wind  profile  with  Lettau's  vertical  mixing  velocity. 

Starting  with  Deacon's  Equation 


^P 


kl  =  u*  /z  +  z»)  (15) 

2     R2a  t   Z„  / 


where  /3  ,  as  defined  by  Deacon,  is  a  decreasing  function  of  the  Richardson 

number.   Combining  Eq.   (7)  and  (8),  one  obtains 
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u.  ,il 


where  u  ■  X   ^ 
3Z 


*     0     *. 
ind  w  =  _^L  vA 


So  that        u 


Upon  substituting  into  Eq.  (15) 


7z      /2+2tf 


<>z  )?zfl 


But    from  Eq .    (14)     ^L   =  * 

4.  0+x) 

the  radical  by  >fe/fe. 


,    so   that   after  multiplying   inside 


du 


.4 


J*,       Wa. 
O+X)    (l+^T 


Z+2C 


r/5 


frz, 


Dividing  through  by   -i/^u 


ylzT 


Squaring   out 


du    =     44         /ztze 


-2.fi 


But    from  Eq.    (2)       ^    -     |?  (z.  *  Z0j       ,    so   that 


^u    _ 


<}Z 


few,rz+zj  /?+ZcX 

**Za2 


or, 


<3a 

c)Z 


fcz, 


,o+x)z[  z0/(   zfly 


-2/9 


Further  combining  of  terms  gives  the  equation  in  its  final  form.   Namely, 


az        fcz^o+xr 


P^j 


/-2/5 


(16) 


Thus  Martin  has  eliminated  u^  from  Deacon's  Equation  -  a  parameter  that 
is  difficult  to  determine  with  accuracy. 
Rewriting  Eq .  (16)  in  the  form 


2 


1-/3 


lez, 


(^) 


_/? 


(17) 


and  comparing  it  with  Eq .  (15),  one  can  easily  see  that  an  equation  for 
the  friction  velocity  has  evolved: 


1-/3 


^(^) 


Using  this  equation  to  find  u.  one  needs  only  to  know  the  adiabatic  mixing 
velocity,  roughness  parameter,  vertical  temperature  gradient,  and  the 
stability  parameter. 

This  completes  the  necessary  derivations.   The  remainder  of  the  paper 
deals  with  the  solution  of  Eq.  (  2)  for  x,  Eq.  (16)  for   /3    ,  and  Eq .  (18) 
for  u^,  using  available  micrometeorological  data. 


3.   Data 

Immediately  upon  deciding  to  conduct  this  investigation  into  micro- 
meteorology  a  search  was  initiated  to  find  suitable  data.   It  soon  became 
apparent  that  there  was  a  scarcity  of  micrometeorological  data  of  the 
caliber  necessary  to  conduct  individual  research  into  surface  layer  turbu- 
lence.  Several  sets  of  available  published  data  had  to  be  discarded 
because  one  or  another  necessary  parameter  had  not  been  observed,  the  period 
of  observation  was  too  short,  or  the  heights  of  the  observations  were  not 
at  suitable  locations. 

The  best  data  that  could  be  found  within  a  reasonable  time  was  that  by 
Gerhardt  [5]. in  1949  and  by  Gerhardt  [6]  in  1950.   The  1949  data  that  was 
of  interest  in  this  investigation  consisted  of  wind  observations  at  12', 
41',  and  91';  aerated  temperatures  at  3',  6',  12',  20',  35',  55',  80',  and 
110';  non-aerated  temperatures  at  1",  6",  18",  3',  5',  9',  and  15'.   The 
1950  data  included  in  addition  to  the  above,  wind  at  3'  and  non-aerated 
temperatures  at  23',  33',  45',  60',  78',  and  99 ' . 

The  information  was  gathered  during  runs  of  from  24  to  44  hours  for 
different  seasons  of  the  year,  with  an  average  run  being  about  30  hours. 
Observations  were  taken  every  half  hour  during  a  run,  with  the  final  value 
for  that  time  and  height  being  an  average  of  five  readings  taken  over  a 
six-minute  interval. 

Jehn  [9]  and  Gerhardt  [7]  give  a  detailed  explanation  of  the  instruments 
used  and  their  accuracy,  the  topography  of  the  micrometeorological  site 
area,  and  the  vegetation  at  the  site  and  surrounding  fields.   The  wind 
measurements  were  taken  with  Bendix-Friez  Aerovanes  to  the  nearest  mph. 
The  errors  involved  were  less  than  1  mph  for  wind  speeds  in  the  range 
10-50  mph,  increasing  to  nearly  2  mph  after  13  months  of  operation.   For 
wind  speeds  less  than  10  mph  the  errors  were  slightly  greater,  and  readings 
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below  3  mph  were  unsatisfactory. 

The  wind  measurements,  because  of  their  inaccuracy  and  the  spacing 
of  instruments,  were  a  severe  handicap  in  this  investigation.   However, 
it  was  hoped  that  by  analyzing  a  large  number  of  observations,  the  errors 
would  be  somewhat  minimized. 

The  temperatures  were  obtained  from  calibrated  rod-shaped  thermistors 
(ceramic  resistors),  about  0.75  inches  long  and  0.03  inches  in  diameter. 
All  thermistors  were  fitted  to  a  single  calibration  curve  giving  an 
accuracy  over  the  range  of  0  C  to  50  C  of  -  0. 1  C.   The  aerated  thermistors 
were  mounted  in  anti-radiation  housings,  being  ventilated  by  a  conven- 
tional motor  and  exhaust  fan  arrangement  at  the  top  of  the  housing.   The 
non-aerated  thermistors  were  coated  with  an  anti-radiation  paint  and,  in 
later  runs,  were  shielded  from  both  solar  and  terrestrial  radiation  without 
reducing  the  natural  ventilation  available. 

Because  of  all  the  precautions  taken  in  calibration  and  shielding, 
•ne  could  expect  that  there  would  be  little  or  no  difference  between  the 
aerated  and  non -aerated  temperatures.   However,  during  this  investigation 
it  soon  became  evident  that  there  were  large  differences  (in  some  cases 
nearly  1  C)  between  the  two  elements,  and  in  several  instances  one  gave  an 
increase  in  temperature  with  height  while  the  other  gave  a  decrease  over 
the  same  layer.   In  order  to  determine  which  would  be  best  suited  for  this 
Investigation,  a  thorough  analysis  of  both  sets  of  temperatures  was  con- 
ducted for  two  runs.   It  was  found  that  at  elevations  greater  than  five 
feet  the  non-aerated  temperatures  were  higher  than  the  aerated  both  day  and 
night,  with  the  differences  usually  .2-. 6  C.for  all  pairs  of  thermistors. 

For  elevations  lower  than  five  feet,  however,  it  was  discovered  that  the 
non-aerated  thermistors  were  cooler  during  the  early  afternoon  and  warmer 

during  the  early  morning,  with  the  differences  increasing  with  decreasing 
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elevation.  This  at  first  seemed  puzzling,  since  it  is  contrary  to  what 
one  vould  intuitively  expect  --  at  least  for  early  afternoon.  However, 
the  author  believes  that  this  can  be  reasonably  explained. 

The  magnitude  of  the  vertical  temperature  gradient  decreases  rapidly 
with  height  near  the  ground,  especially  during  early  morning  and  afternoon. 
Now  consider  the  ventilation  system  drawing  in  an  equal  volume  of  air  from 
both  above  and  below  the  level  z.   Because  of  the  difference  in  tempera- 
ture gradient  the  absolute  value  of  the  temperature  difference  between 
z  and  z  -  4 z  will  be  greater  than  that  between  z  and  z  +  &  2 .   This  dif- 
ference will  be  greatest  during  early  morning  and  early  afternoon,  and 
thus  the  aerated  temperature  will  be  warmer  than  the  true  temperature  at 
z  during  the  afternoon  and  cooler  in  the  morning. 

Because  of  this  effect,  it  was  decided  that  the  non-aerated  thermistors 
would  give  a  more  representative  vertical  gradient.   However,  of  the  four 
runs  used  in  this  paper  (Runs  8,  10,  16,  and  17)  only  on  Run  17  was  it 
possible  to  use  the  non-aerated  temperatures  since  on  Runs  8  and  10  they 
were  only  observed  to  15'  and  on  Run  16  the  3'  element  was  broken.   Even 
so,  it  is  felt  that  greater  accuracy  in  the  temperature  measurements  are 
not  warranted  because  of  the  large  errors  involved  with  the  wind 
velocities. 

Gerhardt,  in  his  final  evaluation  of  the  Micrometeorological  Research 

Project  (from  which  the  data  used  in  this  investigation  was  obtained), 

stated: 

The  importance  of  an  accurate  wind  gradient,  particularly  in  the 
lowest  layers,  cannot  be  overestimated;  12  measurement  levels  below 
SO  feet  would  not  be  excessive.   Great  care  should  also  be  taken  lo 
obtain  Instruments  with  as  nearly  identical  response  as  possible. 
Although  thermistors  have  certain  undesirable  characteristics,  it  is 
believed  that  the  installation  as  set  up  was  satisfactory,  and  that 
accurate  and  representative  air  temperatures  were  obtained. 
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4.   Computations 

In  this  section,  the  working  equations  using  finite  differences  in 
place  of  derivatives  are  illustrated.   The  evaluation  of  the  supposedly 
"non-variable"  parameters  is  discussed,  and  the  advantages  and  dis- 
advantages in  using  an  electronic  digital  computer  for  this  type  of  com- 
putation is  pointed  out. 

The  particular  working  equation  used  to  compute  x  was  obtained  by  sub- 
stituting Eq.  (2)  into  Eq .  (12)  and  substituting  finite  differences  in 
place  of  the  partial  derivative  cXyd2-     .      This  results  in  the  working 
equation 


X  = 


\frir*,) 


62.-6, 


(19) 


The  terms  inside  the  bracket  are  constant  for  the  layer  during  most 
of  the  run,  although  a  new  value  of  w  was  selected  every  12  hours.   It 

0 

is  easily  seen  from  this  equation  that  x  is  zero  for  neutral  conditions, 

greater  than  zero  for  unstable  conditions,  and  less  than  zero  for  stable 

2 
stratification.   Values  of  (1+x)   were  tabulated. 

The  next  variable  to  find  was  /9   .   The  partial  derivative  <)u/^Z  was 

first  calculated  from  the  wind  data  using  finite  differences.  Then  taking 

logarithms  of  both  sides  of  Eq.  (16)  and  solving  for  /?     resulted  in  the 

following  equation. 


/*-  ^ 


z+z* 


1+T-o 


(20) 
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The  only  variable  for  a  given  level  Is  the  numerator  of  the  third  term 
inside  the  brackets. 

The  final  variable  to  determine  for  this  investigation  was  the 
friction  velocity.   The  working  equation  was  developed  by  taking  logarithms 
of  both  sides  of  Eq.  (18)  and  then  taking  e  to  the  power  of  the  left  side 
and  the  right  side,  giving 


.         aA_   +    o-/3) a  ***a  (21) 


This  equation  looks  more  difficult  to  solve  than  the  original  equation, 
but  is  actually  more  adaptable  to  programming  on  our  electronic  digital 
computer. 

Prior  to  solving  Eqs.  (19),  (20),  and  (21)  the  "constant"  parameters 
(von  (Carman's  constant  k,  adiabatic  mixing  velocity  w  ,  and  roughness 
parameter  z  )  had  to  be  determined. 

According  to  [8],  w  and  z   are  easily  obtained  by  plotting  the  wind 
speeds  against  In  z  for  adiabatic  conditions,  and  drawing  a  best-fitting 
straight  line.   The  roughness  parameter  being  the  intercept  on  the  In  z 
axis  and  the  mixing  velocity  being  the  slope  divided  into  von  Karman's 
constant,  or 

2"/ 

The  first  problem  to  be  solved  was  that  of  determining  the  time,  and 
therefore  the  wind-sounding,  for  neutral  conditions.   All  of  the  layers 
examined  up  to  110  feet  were  of  depth  not  exceeding  10  meters.   In  layers 
of  such  depth,  bounded  by  thermistors,  a  dry  adiabatic  lapse  rate 
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corresponds  to  a  vertical  temperature  difference   £~T|  —0.1  C.   In  such 
shallow  layers,  a  neutral  lapse  may  be  regarded  as  essentially  iso- 
thermal.  Special  emphasis  in  this  regard  was  given  to  the  temperature  dif- 
ference in  the  layer  0-15  feet  above  the  ground.   Runs  16  and  17  corre- 
spond to  stabilities  nearly  neutral  for  the  entire  period. 

On  clear  days,  a  neutral  surface  layer  is  considered  [8]  to  occur 
near  sunrise  or  sunset.   This  was  found  to  be  Che  case.   The  actual  time 
of  lsothermalcy  usually  occurred  between  half-hourly  wind  observations  near 
sunrise.   Hence  determining  appropriate  wind  values  required  the  averaging 
of  at  least  two  consecutive  sets  of  data  centered  about  the  time  of  lso- 
thermalcy.  Actually  it  was  decided  to  average  wind  speeds  for  an  hour 
before  and  after  the  neutral  case  in  order  to  reduce  round  off  errors  in 
the  reported  wind  values,  since  these  were  given  only  to  the  nearest  mph. 
These  averaged  wind  speeds  were  then  plotted  against  In  z  and  a  straight 
line  fitted.   The  straight  line  went  through  or  nearly  through  all  the 

plotted  points.   The  values  of  z  were  consistent  with  those  12  hours 
r       r  o 

earlier  or  later.   The  adiabatic  mixing  velocity  w  still  varied  from  the 

a 

morning  to  evening  adiabatic  case.   This  is,  however,  to  be  expected.   The 
adiabatic  mixing  velocity  thus  obtained  was  used  for  six  hours  preceding 
and  six  hours  after  the  time  of  neutral  stratification. 

Suggestions  have  been  made  [16],  [17]  that  k  is  a  function  of  stability 
However,  most  specialists  in  the  field  treat  k  as  a  constant  independent  of 
stability.   Values  given  for  it  in  textbooks  and  in  the  literature  range 
from  0.38  to  0.45,  with  most  references  giving  k  as  0.40.   This  is  the 

value  used  in  this  thesis.   The  particular  constant  value  used  for  k  will 

2 
not  materially  affect  the  results  of  correlating  o   with  (1+x)  . 

# 
It  was  decided  to  use  the  CRA  102-A  general  purpose  automatic  digital 

computer  [19]  rather  than  a  slide  rule  or  hand  computer.   This  decision  was 
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based  on  the  large  number  of  computations  involved  and  the  additional 
numbers  of  significant  digits  possible  with  the  computer. 

There  are  two  possible  disadvantages  in  using  such  a  computer,  the 
first  being  that  the  data  had  to  be  converted  from  the  decimal  number 
system  to  the  octal  number  system.   If  the  data  were  whole  numbers,  such 
as  the  wind  observations,  this  presented  little  difficulty  because  programs 
are  available  for  the  necessary  conversion.   However,  when  involved  with 
mixed  numbers,  as  in  the  temperature  observations,  it  is  necessary  to  con- 
vert the  whole  number  part  and  the  fractional  part  separately,  and  then 
add  them  together.   It  was  possible,  however,  after  a  little  experience  to 
rapidly  convert  the  whole  numbers  mentally,  thus  necessitating  machine 
conversion  of  the  fractions  only.   It  was  found  that  for  a  run  of  24  hours 
with  wind  and  temperature  at  four  heights  and  the  mean  temperatures  of  the 
three  layers  (a  total  of  528  numbers)  that  the  data  could  be  converted  and 
put  on  IBM  cards  in  less  than  one  hour. 

The  second  disadvantage  was  the  writing  of  the  programs  for  the  com- 
puter.  It  is  estimated  that  a  total  of  20  hours  was  spent  writing  the 

2 
necessary  programs  (Ri,  3u,  (1+x)  ,  ft  ,  u^,  and  a  correlation  coefficient 

d  z 
program),  and  another  20  hours  spent  testing  these  programs  on  the  com- 
puter to  make  certain  they  were  working  properly. 

Once  this  preliminary  work  was  accomplished  the  remaining  work  was 
relatively  simple,  since  the  programs  were  good  for  any  run  once  a  few  con- 
stants were  changed.   An  example  of  the  time  saved  will  be  illustrated. 
Five  friction  velocities  were  computed  by  hand  using  logarithms  and  a  slide 
rule,  taking  nearly  one  hour.   Using  the  computer,  the  friction  velocities 
for  one  24-hour  run  and  three  layers  (244  values)  were  computed,  typed  out 
in  both  octal  and  decimal,  and  IBM  cards  punched  in  15  minutes]   Assuming 
that  the  time  for  the  hand  calculations  could  be  reduced  by  half,  it  would 
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still  take  over  24  hours  or  about  100  times  as  long  as  the  machine. 

By  having  the  computer  punch  IBM  cards  for  all  variables  it  was 
then  an  easy  matter  to  determine  correlation  coefficients  between  the 
various  parameters.   It  took  the  machine  36  seconds  to  compute  the 
correlation  coefficient  between  two  parameters  consisting  of  48  observa- 
tions each. 

The  computer  can  also  plot  the  information  taking  a  little  over  a 
minute  to  plot  48  values  of  two  variables.   This  was  especially  useful 
because  the  relatively  small  range  the  various  parameters  had  would  make 
it  very  tedious  to  plot  by  hand. 

Once  a  program  has  been  successfully  tested,  there  is  little  chance 
of  erroneous  values.   With  hand  calculations,  arithmetical  errors  may 
occur  at  any  time. 

One  last  advantage  that  seems  worthy  of  mentioning  is  the  fact  that 
with  all  the  information  on  cards,  it  would  be  easy  for  another  person  to 
use  in  an  investigation  along  similar  lines  -  such  as  in  determining  the 
heat  flux  or  the  coefficient  of  eddy  viscosity. 
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5.   Results 

A  discussion  of  the  results  of  this  investigation  and  a  comparison 
of  these  results  with  those  of  other  investigators  into  surface  layer 

turbulence  will  now  be  made.   The  numerical  values  of  the  parameters 

2 
(1+x)  ,  /? ,  and  u^  for  each  run  may  be  found  in  tables  in  Appendices  I, 

II,  III,  and  IV.   It  will  be  seen  in  these  tables  that  for  Runs  16  and  17 

three  layers  were  evaluated  -  3  to  12  feet,  12  to  41  feet,  and  41  to 

91  feet.   Although  Deacon's  wind  profile  is  usually  considered  applicable 

only  to  about  41  feet,  a  higher  layer  was  analysed  in  order  to  determine 

the  inconsistencies,  if  any,  from  those  at  lower  elevations.   Runs  8  and 

10  were  used,  even  though  it  was  only  possible  to  get  values  for  the 

middle  layer  (due  to  insufficient  data) ,  because  of  the  much  greater  ranges 

of  stability  than  for  Runs  16  and  17. 

2 
The  diurnal  variation  of  (1+x)  ,  /3  ,  and  u^,  as  well  as  their 

correlation  rfith  one  another  i«  illustrated  in  Fig.  1,  where  the  three 

parameters  for  Run  8  (12  to  41  feet)  are  plotted  against  time.   Run  8  was 

taken  on  27-28  September  with  clear  skies  and  light  to  moderate  winds 

2 
prevailing  for  most  of  the  run.   It  is  seen  that  (1+x)   generally  follows 

the  diurnal  march  of  temperature,  becoming  greater  than  one  after  sunset 

and  less  than  one  after  sunrise,  whereas  G   has  the  opposite  trend.   For 

this  particular  run  [3   was  less  than  one  for  only  a  few  hours  just  prior  to 

midnight,  indicating  the  layer  was  unstable  for  the  majority  of  the  run. 

2 
The  inverse  relationship  between  3  and  (l4x)   can  easily  be  recognized 

in  Figs.  1,  2,  and  3.   The  curves  of  Figs.  1  and  6  have  been  smoothed 

somewhat  using  a  smoothing  equation 
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For  a  closer  examination  of  the  short  period  fluctuations  of  (3  ,    see  the 

2 
Appendix.   A  negative  correlation  between  /?  and  (1+x)   is  to  be  expected 

since  x  is  directly  related  to  the  Richardson  number 

/duf 

m 


2 
The  quantities  (3  and  (1+x)   in  Figs.  2  and  3  are  weighted  means  for   the 

entire  layer  from  3  feet  to  91  feet  for  each  observation  time  of  Run  17. 

The  equation  used  to  determine  this  mean  being 


/3  -    __L      9/?      +  2<?/3      +    50/3 
I  go         13-12  /iz-h  Ui-n 


2  

and  a  similar  equation  for  (1+x)  .   Fig.  3  is  a  scatter  diagram  of  (3    and 

1 

(1+x)  ,  with  the  line  fitted  by  eye.   The  scatter  of  the  points  is  prob- 
ably mainly  the  result  of  the  inaccuracy  of  the  original  wind  measurements 
which  were  given  only  to  the  nearest  mph.   For  Run  17  the  linear 


—     2 

correlation  coefficient  between  /?  and  (1+x)   of  -.54  was  determined.   This 

significant  negative  correlation  is  considered  to  indicate  the  soundness 

of  the  Lettau-Martin  theory. 

In  1948,  Sutton  [18]  deduced  that  the  Richardson  number  must  be 

independent  of  height.   However,  Deacon,  in  the  same  year,  found  from  his 

observations  and  calculations  that  Ri  varied  almost  linearly  with  height. 

Fig.  4  shows  (1+x)   for  each  layer  plotted  against  time  for  Run  16.   A 

2 
decrease  of  (1+x)   with  height  (and  thus  a  decrease  of  Ri  with  height)  is 

2 
very  apparent  in  this  figure.   The  average  value  of  (1+x)   was  determined 

for  each  layer,  using  the  usual  equation  to  find  a  mean 
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Figure  5.       The  linear  variation  of   (l+x)|      with  height  for  Run  l6. 
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The  average  values  were  then  plotted  against  height,  resulting  in  Fig.  5. 

2 
The  linear  decrease  of  (1+x)    with  elevation,  for  this  run,  substan- 

avg 

tiates  Deacon's  conclusions. 

In  comparing  the  values  of  /3  obtained  in  this  investigation  with  those 
of  previous  investigators,  they  were  found  to  fall  within  the  same  range. 
Calder  [l]  had  values  of  fi   from  .75  to  1.4,  Deacon  (1948)  found  that/? 
varied  from  .75  to  1.2,  and  Priestly  and  Sheppard  [15]  determined  the 
range  of  Q   to  be  from  .70  to  1.3.   A  minimum  value  of  /3   equal  to  .84 
occurred  in  this  investigation  during  Run  10,  with  the  maximum  value  of 
1.4  appearing  during  Run  8. 


Investigator 

Min/3 

Max  £ 

Calder 

.75 

1.4 

Deacon 

.75 

1.2 

Priestly  &  Sheppard 

.70 

1.3 

Johnston 

.84 

1.4 

Table  of  ranges  of  /3   as  determined  by  the 
indicated  investigators. 

Table  1 


fi   was  found  to  vary  considerably  from  one  observation  time  to  the 
next,  with  a  maximum  variation  of  around  .2  not  uncommon  for  Runs  8  and 
10.   Longley  [11]  derived  an  equation,  for  /?  that  made  the  usual  assumption 
of  z   constant  with  changes  in  stability  unnecessary.   In  this  equation 

p    depends  only  on  wind  speeds  at  three  levels  if  the  heights  are  related 

2 
by  z   x  z2r  »   z^r  t    where  r  is  any  ratio.   His  equation  was 

24 


1  ^n  r* 


Using  this  equation,  he  computed  values  of  (3   from  Deacon's  data  and  found 
an  average  difference  from  Deacon's  /S  of  from  0.1  to  0.2  depending  on  the 
heights  selected.   He  also  determined  (3  from  other  data  for  various  times 
during  both  day  and  night.   He  found  that  for  daylight  hours  the  range  of 
four  values  of  /3   taken  during  the  same  hour  might  lay  between  0.1  and  0.2, 
and  for  inversion  conditions  the  ranges  may  even  be  greater.   His  final 
conclusions  were  (1)  a  possible  error  of  0.1  or  0.2  does  not  permit  dif- 
ferent types  of  stability  to  be  clearly  distinguished  in  many  cases,  and 
(2)  it  is  impossible  to  calculate  /?  with  sufficient  accuracy  to  warrant 
its  use  in  any  practical  problems.   Deacon,  in  answering  Longley's  article, 
agreed  that  j2   is  difficult  to  measure  accurately,  and  said, 

I  have  never  been  so  optimistic  as  to  suppose  that  in  practical 
problems  it  would  be  profitable  to  attempt  the  routine  measure- 
ments of  /3   in  the  field. 

It  should  be  remembered  that  Longley's  method  of  measurement  was 
based  essentially  on  curve  fitting  and  had  no  way  of  separating  out  the 
effects  of  instability.   Longley's  conclusions  may  no  longer  be   tenable 


2~ 

since  in  the  present  paper,  use  of  (1+x)   accounts  for  a  considerable  part 

of  the  variability  of  j[3  and  of  u^. 

In  perusing  the  available  literature  on  surface  layer  turbulence,  etc., 

no  investigations  into  the  possible  variation  of  Q   with  height  were  found. 

f3   is  generally  considered  constant  with  height,  at  least  in  the  lower 

40  feet.   However,  Deacon  (1948)  felt  there  was  a  possibility  that  {9 

might  vary  slightly  with  elevation,  and  that  if  it  did,  it  should  slightly 

increase.   Fig.  6  indicates  a  possible  variation  of  /?  with  height. 
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Figure  7.   The  variation  of  8  with  the  logarithm  of  height  for  Run  16. 
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/3   ,  from  Run  16,  vas  averaged  for  each  layer  using  the  same  equation 

2 
as  that  for  (1+x)    ,  with  the  results  as  below, 
avg 


/£..  0.9702 

/?  =  1.0040 

/  Us- 

/3  =  1.0271 

It  was  found  that  in  plotting  these  values  against  In  z,    as  in  Fig.  7, 
the  points  could  be  connected  by  a  straight  line.   This  suggests  that  (3 
increases  with  the  logarithm  of  height  in  the  friction  layer.   The  average 
(2   for  each  layer  was  also  determined  for  Run  17.   Here  it  was  found  that 
the  average  j3   for  the  upper  layer  decreased  slightly  from  that  for  the 
middle  layer.   However,  a  straight  line  connecting  /3       for  the  middle 
and  lower  layer  gave  a  similar  slope  as  that  found  for  Run  16.   This  could 
possibly  illustrate  the  fact  that  Deacon's  equation  was  not  applicable  to 
as  great  a  height  for  Run  17  as  for  Run  16,  due  to  a  difference  in  thick- 
ness of  the  surface  layer,  or  it  could  mean  that  the  relationship  of 
with  height  found  from  Run  16  was  only  accidental. 

Prior  to  this  paper,  there  has  been  no  equation  for  accurately  determin- 
ing the  friction  velocity  for  non-neutral  conditions  without  measurements  of 
the  surface  drag  or  short  period  wind  velocity  fluctuations.   Panofsky  [14], 
by  making  several  assumptions,  derived  an  equation  for  u^  which  depends 
only  on  the  wind  speed  at  three  elevations.   However,  upon  comparing  his 
results  for  u..  */ith  those  obtained  by  measurements  of  the  short  period 
velocity  fluctuations,  they  were  found  to  be  less  than  satisfactory. 

Deacon  [3],  while  making  an  investigation  into  the  variation  of  the 
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eddy  stress  with  height,  found  that  values  of  u  +  at  75  feet  were  much 
more  variable  than  those  at  five  feet.   He  attributed  this  mainly  to 
accelerations  associated  with  large  scale  eddies.   He  also  found  that 
there  seemed  to  be  no  systematic  variation  of  u^  with  height,  although 
individual  measurements  showed  large  departures. 

Fig.  8  is  the  friction  velocity  from  Run  16  for  the  three  layers  plotted 
against  time,  as  determined  by  Eq .  (18).   Here  it  is  seen  that  the  fluctua- 
tions of  u^,  for  the  highest  layer  are  generally  of  greater  amplitude  than 
for  the  lower  two  layers.   Although  this  run  seems  to  indicate  a  possible 
increase  of  u^  with  height,  Run  17  indicated  a  possible  decrease  with  height. 
However,  there  does  not  seem  to  be  a  systematic  increase  of  u.  between  the 
layer  3-12  feet  as  compared  to  12-41  feet.   Any  possible  conclusions  that 
can  be  made  here  regarding  the  variation  of  u^  in  tue  layer  41-91  feet  are 
tentative,  inasmuch  as  the  Deacon  profile  may  not  be  applicable  to  these 
elevations.   However  there  is  no  clear-cut  increase  of  u^.  between  the  layer 
12-41  feet  and  41-91  feet. 
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6.   Summary  and  Conclusions 

Martin's  combination  of  Lettau's  vertical  nixing  velocity  with 
Deacon's  wind  profile  equation  gives  values  of  ft   and  u   which  compare 
favorably  with  those  of  other  investigators.   However,  the  equations  con- 
taining  the  term  (1+x)   will  not  apply  for  x<  -1,  which  may  have  the 
significance  of  placing  an  upper  limit  on  the  depth  of  the  surface  layer. 
During  this  investigation  x  was  never  found  to  be  less  than  minus  one. 

The  equation  for  the  friction  velocity,  Eq .  (18),  may  become  an 
accurate  method  for  determining  u+  if  it  becomes  possible  to  evaluate  (3 
with  sufficient  accuracy.   Additional  work  is  needed  to  compare  u^,  as 
determined  from  Eq .  (18),  with  that  obtained  from  surface  drag  measurements 
or  short  period  wind  velocity  fluctuations,  before  an  accurate  evaluation 
of  the  equation  can  be  made. 

By  comparing  a  large  number  of  cases  it  was  found  possible  to  use  data 
that  was  otherwise  undesirable  because  of  the  inaccuracies  of  the  wind 
measurements.   However,  to  determine  the  various  parameters  for  any  one 
observation,  it  is  felt  that  wind  measurements  of  much  greater  accuracy  are 
needed. 

One  of  the  most  important  results  of  this  investigation  is  the  determin- 
ing of  a  possible  variation  of  /?  ->ith  height  and  the  correlation  of  {Z   with 

2 
(1+x)  .   However,  further  work  with  a  larger  number  of  layers  is  needed  to 

verify  this  result. 
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Appendix    I 

Run    8 

time 

(l+x? 

-A- 

U*    Ctr>ph) 

1300 

.528 

1.133 

.899 

1330 

.496 

1.226 

.583 

1400 

.496 

1.226 

.583 

1430 

.497 

1.141 

.919 

1500 

.598 

1.168 

.661 

1530 

.633 

1.161 

.648 

1600 

.708 

k.147 

.623 

1630 

.708 

1.119 

.721 

1700 

.827 

1.163 

.489 

1730 

.911 

1.201 

.362 

1800 

1.000 

1.275 

.223 

1830 

1.399 

1.098 

.408 

1900 

1.684 

1.076 

.383 

1930 

2.196 

1.007 

.422 

2000 

2.336 

.972 

.478 

2030 

'    2.269 

.976 

.483 

2100 

2.487 

.964 

.468 

2130 

2.208 

1.007 

.421 

2200 

1.879 

.976 

.481 

2230 

1.945 

1.022 

.440 

2300 

2.078 

1.014 

.430 

2330 

1.758 

1.007 

.527 

0000 

1.408 

1.062 

.492 

0030 

1.638 

1.043 

.467 

0100 

1.464 

1.030 

.562 

0130 

1.464 

1.093 

.402 
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Run  8  (cont.) 

time 

O+xf 

-Jl- 

u* 

0200 

1.408 

1.098 

.407 

0230 

1.522 

1.088 

.396 

0300 

1.409 

1.097 

.407 

0330 

1.145 

1.173 

.335 

0400 

1.196 

1.118 

.431 

0430 

1.196 

1.118 

.431 

0500 

1.249 

1.077 

.512 

0530 

1.302 

1.157 

.320 

0600 

1.146 

1.123 

.437 

0630 

1.125 

1.067 

.373 

0700 

.882 

1.147 

.311 

0730 

.670 

1.275 

.257 

0800 

.409 

1.327 

.257 

0830 

.492 

1.303 

.178 

0900 

.216 

1.405 

.320 

0930 

.277 

1.375 

.294 

1000 

.497 

1.217 

.379 

1030 

.221 

1.403 

.318 

1100 

.347 

1.212 

.560 

1130 

.223 

1.266 

.652 

1200 

.224 

1.266 

.652 

1230 

.089 

1.379 

.893 
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Appendix   II 

Run    10 

time 

0+xf 

-£- 

^*    (rnph) 

1630 

.459 

1.103 

.627 

1700 

.575 

1.081 

.565 

1730 

.574 

1.081 

.566 

1800 

2.716 

.873 

.368 

1830 

1.112 

.956 

.575 

1900 

2.217 

.854 

.499 

1930 

2.470 

.844 

.472 

2000 

1.552 

.989 

.344 

2030 

.032 

1.150 

1.610 

2100 

.180 

1.124 

1.430 

2130 

2.233 

.891 

.406 

2200 

2.488 

.881 

.384 

2230 

1.557 

.989 

.343 

2300 

2.493 

.945 

.271 

2330 

1.359 

.937 

.520 

0000 

1.173 

.951 

.560 

0030 

1.000 

.966 

.606 

0100 

1.361 

.937 

.520 

0130 

1.174 

1.015 

.396 

0200 

1.440 

.918 

.704 

0230 

1.563 

.910 

.675 

0300 

1.441 

.918 

.704 

0330 

1.324 

.963 

.599 

0400 

1.566 

.948 

.551 
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Run    10    (contd.) 


time 

0+xf. 

1.324 

-A 

.963 

u* 

0430 

.599 

0500 

1.000 

.952 

.845 

0530 

.559 

1.043 

.922 

0600 

.486 

1.056 

.989 

0630 

1.444 

.918 

.703 

0700 

1.325 

.926 

.734 

0730 

.722 

1.019 

.812 

0800 

.423 

1.133 

.750 

0830 

.423 

1.069 

1.060 

0900 

.492 

1.119 

.695 

0930 

.642 

1.095 

.608 

1000 

.165 

1.220 

1.200 

1030 

.257 

1.201 

.856 

1100 

.368 

1.082 

1.137 

1130 

.498 

1.054 

.977 

1200 

.048 

1.089 

4.596 
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